Integral formulation for a circular cylindrical cavity in infinite solid and a finite length coaxial cylindrical crack compressed axially

A method for solving problems of fracture of an infinite solid with a circular cylindrical cavity and a coaxial cylindrical crack near the surface under an uniform axial compression is proposed using a non-classical criterial approach associated with a mechanism of a local stability loss near the defect. The theory of integral Fourier transforms and series expansions are used to reduce these problems to a system of paired integral equations and then to a system of linear algebraic equations with respect to the contraction parameter.     

转向剂形成人工应力遮挡的实验研究及有限元分析

针对水力压裂过程中采用的人工应力遮挡控制缝高技术,通过室内物理模拟实验,在转向剂浓度一定的条件下考虑携带液与地层压差、渗流量、时间等因素,采用电测法对人工岩芯表面应变进行了测试;利用非线性有限元方法对裂缝尖端应力进行计算,研究了转向剂形成人工应力遮挡的作用效果。结果表明采用转向剂可以建立3~4MPa的人工应力遮档,对采用转向剂的压裂设计和施工参数提出了建议。     

From crack deflection to lattice vibrations—macro to atomistic examination of dynamic cleavage fracture

A set of cleavage experiments with strip-shaped single-crystal silicon specimens subjected to three-point bending is reported. The experiments enabled examination of the relationships between the dynamic energy release rate, the velocity, the orientation-dependent cleavage energy, and the cleavage plane of propagation.Dynamic crack propagation experiments show that when a [0 0 1] silicon single crystal is fractured under three-point bending at ‘parallel’ velocity (directly measured at the bottom surface of the specimen) of up to , it prefers to cleave along the vertical (1 1 0) plane, while when the specimen is fractured under the same conditions but at a velocity higher than , it cleaves along the inclined (1 1 1) plane. At intermediate velocities, the crack will deflect from the (1 1 0) plane to the (1 1 1) plane. Crack velocity was determined by the initial notch length. The local (calculated) velocity of deflection between the cleavage planes ranges from , for a crack propagating on the (1 1 0) plane in the direction, to about , for a crack on the (1 1 0) plane, but in the [0 0 1] direction.It is suggested that the cause of the deflection phenomenon is the anisotropic, velocity-dependent cleavage energy, resulted phonon radiation caused by anisotropic, velocity-dependent lattice vibrations. We have studied the effect of material properties and propose selection criteria to explain the deflection phenomenon: the crack will deflect to the plane of least-energy, for which G−Γ_i(V)=max, or to the plane with maximum crack tip velocity, V_i(Γ)=max.     

临界空穴扩张比准则的适用性

用新研制成的DE36平台钢制的轴对称和平面应变试样,进一步考察了临界空穴扩张比V_(Gi)准则[1]、[2]、[3]的适用性。实验和大应变有限元理论计算结果表明,V_(Gi)准则适用于所研究材料。文中指出,按V_(Gi)准则可根据试样或构件材料中V_G的发展变化来预起裂位置。     

Ductile fracture criteria for simulating shear by node separation method

This study is concerned with the mechanics of ductile fracture in bulk metal forming processes by a finite-element analysis and experiments. Developed is a computer program using conventional finite-element method such that the behavior of crack propagation after ductile fracture can be analyzed. The phenomenon of a material separating into two pieces upon shearing and tensile tests has been simulated using the developed computer program. Special attention is focused on the effect of various ductile fracture criteria on crack initiation and propagation during shearing and tensile tests.     

各向异性两相材料尖劈奇性场的非协调元分析

提出了一个基于位移的、分析柱状各向异性两相材料尖劈端部邻域的奇性位移场和应力场问题的非协调元特征分析法. 该方法从柱状扇区的散度定理出发，将柱状扇区控制方程的弱式化为一个与虚功原理相同形式的方程，采用一种新的非协调元技术把所导出的``虚功原理'转化为标准一阶特征方程的求解问题. 非协调元法中，尖劈端部邻域的位移场假定没有采用奇异变换技术，有限元的单元形式是一维的. 将柱状各向异性两相材料尖劈视为``广义平面应变'问题，位移场与坐标z无关，只关注界面端的幂奇异性而不考虑对数奇异性. 运用该方法给出了柱状各向异性两相材料尖劈端部奇性应力指数、奇性位移角分布和应力角分布的算例. 所有的计算结果表明，该方法使用的单元少而且精度较高.     

薄板试样斜断口形成机制及损伤断裂过程分析

为分析不同材料和尺寸的薄板试样在室温下拉伸破坏后均形成与横截面夹角在20°~25°之间斜断口的原因,首先用统计方法对试样内随机分布微缺陷进行讨论,提出一种在宏观尺度上材料内微缺陷分布局部非均匀简化模型的假设.应用含孔材料损伤本构模型对含有不同方向微缺陷分布局部非均匀薄带区域的16MnNb薄板试样变形至破坏全过程进行数值模拟.结果表明,斜断口形成主要是由于试样内在与横截面夹角小于45°的带形区域内微缺陷分布局部非均匀造成,且与该带形区域在试样中位置无关;由于考虑微缺陷分布局部非均匀,得到试样的斜断口形成过程与试验现象完全一致;同时结合试验断口形貌,对变形过程中颈缩截面内损伤演化和破坏过程进行研究,进一步解释薄板试样的损伤破坏机制.     

Mode III fracture analysis by Trefftz boundary element method

This paper presents a hybrid Trefftz (HT) boundary element method (BEM) by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor (SIF) is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by using the Galerkin and the collocation techniques is effective for mode III fracture problems. The project supported by the National Natural Science Foundation of China(10472082). The English text was polished by Keren Wang.     

A material force method for inelastic fracture mechanics

A material force method is proposed for evaluating the energy release rate and work rate of dissipation for fracture in inelastic materials. The inelastic material response is characterized by an internal variable model with an explicitly defined free energy density and dissipation potential. Expressions for the global material and dissipation forces are obtained from a global balance of energy-momentum that incorporates dissipation from inelastic material behavior. It is shown that in the special case of steady-state growth, the global dissipation force equals the work rate of dissipation, and the global material force and J-integral methods are equivalent. For implementation in finite element computations, an equivalent domain expression of the global material force is developed from the weak form of the energy-momentum balance. The method is applied to model problems of cohesive fracture in a remote K-field for viscoelasticity and elastoplasticity. The viscoelastic problem is used to compare various element discretizations in combination with different schemes for computing strain gradients. For the elastoplastic problem, the effects of cohesive and bulk properties on the plastic dissipation are examined using calculations of the global dissipation force.